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Abstract 

We present preliminary results from a lepton-tagged fully-inclusive measurement of B ^ Xg^ de- 
cays, where Xg is any strange hadronic state. Results are based on a BABAR data set of 88.5 
million BB pairs at the T(45) resonance. We present a reconstructed photon energy spectrum in 
the T{4:S) frame, and partial branching fractions above minimum reconstructed photon energies of 
1.9, 2.0, 2.1 and 2.2 GeV. We then convert these to measurements of partial branching fractions 
and truncated first and second moments of the true photon energy distribution in the B rest frame, 
above the same minimum photon energy values. The full correlation matrices between the first and 
second moments are included to allow fitting to any parameterized theoretical calculation. We also 
measure the direct CP asymmetry Acp{B — > Xg+d^) (based on the charge of the tagging lepton) 
above a reconstructed photon energy of 2.2 GeV. 
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1 Introduction 



In the Standard Model the radiative decay 6 — > 57 proceeds via a loop diagram, and is sensitive to 
possible new physics, with new heavy particles participating in the loop [Tl. Next-to-leading-order 
calculations for the branching fraction give B{B Xs-f) = (3.57 ± 0.30) x 10"^ {Ej > 1.6 GeV) 
[2], where is the photon energy in the B rest frame; calculations to higher order are currently 
underway [2]. Measurement of the E^ spectrum from B — > Xg^ decays, and in particular of its 
moments above various minimum energies, allows the determination of the heavy quark effective 
field theory parameters nib and related to the 6-quark mass and momentum within the B 
meson, respectively jH El El El- 

This note reports on a fully inclusive analysis of data collected from T{4S) BB, where the 
photon from the decay of one B meson is detected, but the Xg system is not reconstructed. Previous 
fully-inclusive measurements of i? — > X^j have been presented by the CLEO (8j and BELLE ^ 
collaborations. The alternative approach of summing a number of exclusive decays jlOj is more 
dependent on the Xg fragmentation model, and in particular on assumptions made as to the fraction 
of unmeasured final states. In the fully-inclusive measurement we largely avoid such uncertainties, 
but at the cost of higher backgrounds, which need to be strongly suppressed. Much of the non-BB 
(continuum) background is removed by requiring a high-momentum lepton tag; this selects BB 
events in which the non-signal B meson decays semileptonically. Continuum background is further 
suppressed by topological cuts. Remaining continuum background is subtracted using off-resonance 
data recorded at an e~^e~ center-of-mass (cm.) energy just below that of the T(4S'). Backgrounds 
from BB events are estimated from a Monte Carlo simulation; for the most important backgrounds 
the simulation is checked and corrected using data control samples. 

We present preliminary results of this analysis, in particular the partial branching fraction and 
the first and centralized second moments of the photon energy distribution {{E^) and {E^) -{E^f, 
respectively) for each of four minimum- i?^ cuts in the B rest frame: 1.9, 2.0, 2.1 and 2.2 GeV. 
Hereafter we refer to the centralized second moment as the "second moment". These are obtained 
from the corresponding quantities for reconstructed photon energy E* in the e~^e~ cm. frame (an 
asterisk denotes a a cm. quantity) by applying corrections for resolution smearing and the effect of 
imposing an experimental cut on E* rather than E^. The full correlation matrix between the first 
and second moments is given to facilitate fitting to any parameterized theoretical calculation. The 
measurement does not distinguish B X^rf events from B — s- Xg'^ events. We assume the ratio of 
B — > Xrf7 io B ^ Xg'^ is |Vtd/V(c,p, so that it constitutes (4.0 ± 1.6)% of the rate and correct 
for this at the end. This assumption is supported by newer theoretical calculations |111 112j . We 
also present a measurement of the CP asymmetry in the sum ol B ^ Xg^ plus B X^^, based 
on the charge of the tagging lepton. Pending a more detailed study of the energy spectrum, we do 
not yet perform the model-dependent extrapolation to the entire photon spectrum. 

The results presented are based on data recorded with the BABAR detector ^31 at the PEP- 
II asymmetric-energy e~^e~ collider located at the Stanford Linear Accelerator Center. The on- 
resonance data integrated luminosity is 81.5 fb~^, corresponding to 88.5 million BB events. Ad- 
ditionally, 9.6 fb"^ of off- resonance data are used in the continuum background subtraction. The 
BABAR Monte Carlo simulation program, based on GEANT 4 [HI and JETSET [13 , is used to gen- 
erate samples of signal events, B^B^ and B^B^ (excluding signal channels), qq {udsc continuum) 
and T^T~ . In order to minimize any possible experimenter's bias, the spectrum of reconstructed 
photon energy E* from 1.9 to 2.9 GeV in the on-resonance data was "blinded", i.e., not looked at 
until all selection cuts were set and dominant BB corrections determined. 
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2 Signal Model and Backgrounds 



In order to compute signal {B Xgj) efficiency, we use the mx^ spectrum computed by Kagan 
and Neubert (KN) This quantity has a one-to-one correspondence with the photon energy in 
the i?-meson rest frame via the kinematic relation 

2 2 

mi, — niyr 

P = E± . (1 

^ 2mB 

The KN spectrum is a smooth function of mnxs-, J-c, it does not explicitly model strange resonances 
in the mass spectrum. This assumption is not valid for the decay B i^'* (892)7, so the spectrum 
needs to be modified to describe it. Higher- mass resonances are broad and spaced closely enough, 
such that the smooth approximation of the KN treatment is justified. In practice, a relativistic 
Breit-Wigner distribution modeling the B K*^ decay is used for mx^ < 1.1 GeV/c^ (ii^^ > 
2.525 GeV) and the KN spectrum is used above the cutoff. The K*^ part is normalized to the 
fraction of the KN spectrum below the mx,, cutoff. The KN model is parameterized in terms of 
two quantities, (the 6-quark mass) and /x^ (a measure of the quark's Fermi momentum inside 
the B meson). We chose a KN480 model, with mf, = 4.80 GeV/ and fi^ = 0.3 GeV^, for studying 
systematic effects on signal efficiency and for optimizing selection criteria. For the latter purpose, an 
overall normalization was applied such that the B — > K*^ component matches the previous BABAR 
measurement In order to assess possible model-dependence of our results we use additional 
KN models with between 4.55 and 4.80 GeV/c^, along with a set of models based on the theory 
of ref . jSj . 

Backgrounds arise from continuum (e"*"e~ qq, where q = udsc, and r+r^) and BB events. 
There are additional contributions from QED and hadronic two-photon processes, which are not 
simulated. In BB events, the main backgrounds come from decays of vr*^ and rj mesons, with 
a smaller contribution from hadrons (mostly antineutrons) that are misidentified as photons. In 
continuum events, these backgrounds also contribute, along with initial state radiation (ISR), where 
a high-energy photon is radiated from one of the colliding e^. 



3 Event Selection 

The BABAR detector ^S] tracks charged particles with a five-layer double-sided silicon-strip detector 
and a forty-layer drift chamber placed inside a 1.5 T solenoidal magnet. Electromagnetic showers are 
detected in a total-absorbtion calorimeter consisting of 6580 CsI(Tl) scintillating crystals arranged 
hermetically in a projective geometry. Muons are indentified by resistive-plate chambers interleaved 
in the flux-return yoke of the magnet, outside of the calorimeter 

The distinctive feature of the signal decays is an isolated photon with cm. energy from about 
1.6 GeV to 2.8 GeV. Our starting point is finding in the event at least one isolated photon candidate, 
i.e., a localized calorimeter energy deposit, with 1.5 < E* < 3.5 GeV ("photon acceptance"). We 
then apply several categories of selection cuts, described in the next few paragraphs: hadronic 
event selection, photon quality criteria, event shape cuts and tagging requirements. Tabled shows 
the efficiencies of the different categories of cuts for signal and background event samples. In order 
to avoid regions where the expected signal would still be overwhelmed by background (especially 
BB background at low E*), our results will be quoted for several minimum- i?* cuts of 1.9 GeV or 
above. The table illustrates the efficiencies for one particular choice. 

Hadronic events are selected by requiring at least three reconstructed charged tracks and the 
second Fox- Wolfram moment R2 < 0.9, calculated in the cm. frame using all detected particles. 
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To reduce QED and two-photon backgrounds, we use a cut on "effective multiplicity," defined as 
the number of charged tracks plus half the number of photons with energy above 0.08 GeV. We 
require Ues ^ 4.5, a value that has essentially no impact on signal efficiency. 

Photon quality cuts are designed to suppress high-energy photon candidates which arise from 
backgrounds. We require that the polar angle of the photon satisfy —0.74 < cos 9^ < 0.93, to 
ensure that the photon shower is well contained in the calorimeter. We suppress hadronic showers 
and TT^s with merged photon showers by requiring that the shower has a lateral profile consistent 
with a single photon. Photons that are consistent with originating from an identifiable vr'^ or rj 
decay are explicitly vetoed. The optimal TT^{r]) veto requirements are energy-dependent, with a 
more open selection for E* > 2.3 GeV. Finally, isolated photons are selected by requiring that no 
energy (neutral or charged) be within 25 cm (approximately seven times the Moliere radius for 
Csl) of the photon, as measured at on the calorimeter's inside surface. 

Event shape variables are employed to exploit the difference in topology of BB events and 
continuum events. The former tend to be spherical in shape, since the B mesons are produced 
nearly at rest in the T(4S) rest frame, while continuum events tend to have a more jet-like structure. 
We combine 19 different variables into a single Fisher discriminant. Eighteen of the quantities are 
the sum of charged and neutral energy found in 10-degree cones (from to 180 degrees) centered 
on the photon candidate direction. The photon energy is not included in any of the cones. The 
nineteenth variable is the quantity R'^/R^-, where R'2 is the second Fox- Wolfram moment calculated 
in the frame recoiling against the photon (which for ISR events is the qq rest frame), using all 
particles except the photon. The Fisher coefficients were determined by maximizing the separation 
power between simulated signal and continuum events. Cuts on the Fisher variable and i?2 (< 0.55) 
remove 84% of the remaining continuum background at this stage. 

Lepton tagging is used to further reduce the backgrounds from continuum events. About 20% of 
B mesons decay semileptonically to either e or /u, and the large mass of the B meson tends to impart 
large momentum to its lepton daughter. Leptons from fragmentation in continuum events tend to 
be at lower momentum. Since the tag lepton comes from the recoil B meson, it does not compromise 
the inclusiveness of the B Xg'j selection. Electrons are identified by requiring a well-measured 
track matched to an energy deposit in the calorimeter consistent with an electromagnetic shower. 
Muons are identified by matching a well-measured track to hits in the muon detector, with no more 
than a minimum-ionizing energy deposit in the calorimeter. The tag lepton is required to have cm. 
momentum greater than 1.25GeV/c for electrons and 1.5GeV/c for muons. Furthermore, making 
a cut on the angle between the photon and the lepton, cos 9-yi > —0.7, removes more continuum 
background, in which the lepton and photon candidates tend to be back-to-back. Finally, we exploit 
the presence of a relatively high-energy neutrino in semileptonic B decays and place a cut on the 
missing energy in the event in the cm. frame, -E^jsg > 0.8 GeV. 

Virtually all of the tagging leptons arise from the decay B Xctv. We correct the rate of 
such events in our simulation based on the studies presented in Ref. jl7j . The overall rate of 
B — > Xc^v is comprised of four decay modes, where the hadronic system is D, D* , D** (higher 
mass resonances), and non-resonant D^*^7r, respectively. A scale factor was determined for each 
mode and we have weighted our simulated semileptonic decays accordingly. We also reweight the 
simulated form factors of the D*£u decay to those measured by BABAR[TS]. 

The values for many of the cuts described above, along with the Fisher coefficients, were opti- 
mized with all other cuts applied, using an iterative technique. In each case, the quantity maximized 
was S^/(S + B + aC), where S is the expected signal from the KN480 model described above, B 
is the simulated BB background, and C is the simulated continuum background, each being the 
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09'^ 


194 


Hadronic 


U.yo r 


o.y (y 


O.OZD 


Photon quality 


U. / Oo 


U. i ( 


n 1 Qfi 

u. lyo 


2.0 <E*< 2.7 GeV 


0.910 


0.114 


0.326 


Event shape 


0.601 


0.401 


0.167 


Tagging 


0.052 


0.041 


7.4 X 10-"^ 


Total selection 


2.0 X 10-2 


3.2 X 10"^ 


6.5 X 10-6 



Table 1: Selection efficiencies for signal and background events, broken down into different cut cat- 
egories, as described in the text. The total selection efficiency is relative to the photon acceptance. 
The signal model used is KN480 including -fC*7, as defined in the text. The error on the total signal 
efficiency from Monte Carlo statistics is 1.1%. 

number of events expected for the on-resonance data luminosity, and a is the ratio of total to 
off-resonance luminosity. The factor a takes into account the fact that we later use off-resonance 
data to subtract the residual continuum background. 

Figure n shows the overall signal efficiency in bins of E*. The E* dependence is a result of 
the optimization procedure, in particular the training of the Fisher discriminant. The sum of the 
18 energy-flow variables combined in the Fisher discriminant is correlated with the energy of the 
photon in the B rest frame. Without the Fisher selection the continuum background decreases 
by a factor of two over the 2.0 to 2.7 GeV range of E* used for optimization, while the expected 
signal increases rapidly betewen 2.0 and about 2.4 GeV. Therefore, the optimization suppresses 
background more strongly with decreasing i?*, leading to a corresponding reduction of the efficiency. 
The more stringent T^^iff) vetoes applied in events with E* < 2.3 GeV also contribute to the energy 
dependence. The variation of the efficiency with E* does not introduce a significant dependence 
on the multiplicity of the Xg system, which would lead to uncertainties related to fragmentation 
modeling. However, the variation of the efficiency with E* does lead to a dependence on the 
modeling of the photon energy spectrum. The dominant effect is that reducing the assumed value 
of mf, in a particular model leads to reduction in the mean E*, which in turn yields a smaller 
selection efficiency. We improve the estimation of the signal efficiency and reduce the model- 
dependent uncertainty with a bootstrapping procedure described in section 

The simulated E* distribution of signal, continuum and BB background events after all se- 
lection criteria is shown in Fig. |21 Note that the continuum Monte-Carlo does not include QED 
backgrounds present in the data. The actual continuum background is subtracted using the off- 
resonance data and is the dominant source of statistical uncertainty on the measurements. 

4 BB Backgrounds 

As can be seen in Figure [21 the BB background dominates for E* < 2.4 GeV. We present in 
Table [21 the various sources of background coming from BB events, as predicted by our Monte 
Carlo simulation. 

About 80% of the background comes from photons from '7T^{r]) decays, for both the nominal 
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Figure 1: The overall signal efficiency in bins of E* (computed using the KN480 model as defined 
in the text). 

signal region (2.0 < E* < 2.7 GeV) and a control region (1.6 < E* < 1.9 GeV). Other significant 
sources are hadrons in the signal region 8%) and electrons in the control region (~ 9%). We 
have performed dedicated studies of data control samples for these sources, to verify and/or correct 
the Monte Carlo predictions. 

The level of the largest BB background component, T^^irj) decays, has been measured directly 
from the data. We explicitly measure the yield of vr'' and r] decays in bins of E*o or E* in on- 
resonance data, off-resonance data and BB simulation samples, by fitting the 77 mass distributions. 
This allows us to derive correction factors versus tt^{i]) energy for the TT^ii]) rates predicted by the 
simulation. These corrections are then applied to our BB simulation when estimating the BB 
background passing the B — > Xsj selection. The Tr^{i]) analysis uses the same selection criteria 
as for the B Xsj analysis, with three exceptions: 1) the vr''(r/) vetoes are not applied, 2) the 
requirement on the photon energy is reduced to E* > 1.0 GeV to increase statistics, and 3) the 
tagging momentum cuts are loosened to l.OGeV/c for electrons and l.lGeV/c for muons, again to 
increase the available statistics. We apply an additional small correction to account for differences 
in the performance of the vr'^(ry) veto in data and BB simulation. 

The background from hadrons faking photons in BB events consists primarily of ns with a 
small contribution from ps. We have corrected the BB simulation for both the n(p) production 
and n(p) response in the calorimeter. First, an inclusive measurement of p and p production, 
is used to correct the simulation prediction for both p and n production. A conservative 50% 
systematic uncertainty is assigned to account for the fact that we are using a p measurement to 
correct n production. Additionally, A pTT~^ decays are studied both in real data and simulation 
to compare the response of the calorimeter to ps. We find a significantly smaller number passing 
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Figure 2: Monte Carlo simulation of the E* distribution for B Xg^, continuum {uu, dd, ss, cc, 
T~^T~) and BB background after all selection criteria (except that on E* itself) have been applied. 

the E* and lateral-shower-profile cuts in the data than that expected from the simulation. The 
total correction factor, then, is a product of the production correction and the efficiency correction, 
determined in bins of E*. 

Electrons can cause a background to the photon sample if the charged track is not reconstructed 
or matched to the electromagnetic shower. (Bremsstrahlung photons are an additional source of 
background from e^, which is well modeled in the simulation.) The electron tracking efficiency 
was measured in data by selecting one track of a Bhabha event and then measuring the probability 
of reconstructing the other electron track. A further correction is made to account for the fact 
that the environment of the control sample (Bhabha events) has significantly lower multiplicity on 
average than the actual data sample (hadronic events). 

We have also derived corrections for the very small number of background photons from uj and 
ri' decays. The corrections to the modeling oi B —i- Xclf described in section |31 are applied to the 
BB simulation. After all corrections, we check our estimate of the BB background by considering 
the low energy control region 1.6 < E* < 1.9 GeV, where we expect few photons from B — > Xg^ 
decays. In this region we observe 1790 zt 64 events after continuum subtraction compared to an 
expection of 1667 ± 54 events, i.e. 

Observed - Expected = 123 ± 64(stat.) ± 54(syst.) (2) 

The significance of the deviation is l.So". Note that for this check we have not subtracted any signal 
contribution, which for the KN480 and KN465 models defined above would be approximately 17 
and 39 events, respectively. These would reduce the significance of the deviation to l.So" and l.Oo", 
respectively. 
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Fraction of BB background 


Truth-match 


Parentage 


(2.0 <E* < 2.7 GeV) 


Photon 


7\ 


0.640 






0.174 




UJ 


0.024 




r]' 


0.011 




FSR 


0.007 




J/ilj 


0.008 




Other 


0.001 




Total 


0.865 


vr^ (merged) 


Any 


0.001 


Electron 


Any 


0.036 


n [n) 


Any 


0.077 


P {p) 


Any 


0.005 




Any 


0.001 


7r± or 


Any 


0.001 


Unidentified 




0.015 



Table 2: Breakdown of Monte Carlo simulated BB backgrounds by high-energy photon origin for 

2.0 < E* < 2.7 GeV. These fractions are for the combined B^B^ and B^B~ Monte Carlo samples 
before tagging cuts. 

5 Results 

5.1 Partial Branching Fractions 

The extracted signal yield, i.e., the on-resonance yield minus the sum of the scaled off-resonance 
yield and the corrected BB simulation prediction, is given in Table|21for each of four different energy 
ranges. The corresponding raw photon energy spectrum, uncorrected for efficiency, is shown in 
Figure The region 2.9 < E* < 3.4 GeV serves as a control region for the continuum subtraction. 
In this region we observe 390 it 20 events compared to an expectation from off-resonance data of 
391 lb 57 events; a difference of —1 it 61 events. 

The model-dependence of the efficiency can be reduced by a bootstrapping method, whereby we 
initially assume a model (e.g., KN465) and use it to calculate the mean of the efficiency-corrected 
energy spectrum between 2.0 and 2.7 GeV. This is compared to the corresponding mean predicted 
by several KN models with different values of mi,- An iterative procedure is used to determine which 
model best matches the measured mean and to then assign a model-dependent uncertainty to the 
efficiency corresponding to the uncertainty on the measurement of the mean. The relation between 



E; (GeV) 


Tagged Signal Yield (events) 


1.9 to 2.7 

2.0 to 2.7 

2.1 to 2.7 

2.2 to 2.7 


1042 ± 84 ± 62 
992 ± 77 ± 53 
895 ± 72 ± 45 
758 ± 66 ± 40 



Table 3: Extracted tagged signal yields in various ranges of E*, uncorrected for efficiency. The first 
error is statistical, the second is from BB systematics. 
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Figure 3: Photon energy spectrum after background subtraction, uncorrected for efficiency. Error 
bars include statistical (dominant) and BB systematic uncertainties, added in quadrature. 



efficiency and mean energy is approximately linear, with a best match close to a KN460 model 
(with //^ = 0.603). We have carried out the same studies for twelve different parameterizations of 
the Benson, Bigi and Uraltsev (BBU) calculation, using the same prescription to include a K* 
component as for KN. We varied nib between 4.45 and 4.75 GeV and //^ between 0.25 and 0.65 GeV^. 
(Note that these two parameters are defined differently than their KN counterparts.) We set 
(the matrix element of the chromomagnetic operator, responsible for the hyperfine splitting between 
B and B* mesons) to either 0.35 or 0.27 GeV^. All the BBU results for the signal efficiency versus 
mean energy lie on the same line as for the KN model, confirming the model-independence of the 
bootstrapping procedure. Thus the KN460 efficiency is used to determine the partial branching 
fractions (PBFs) according to: 

Kin V TP* ■ \ Events in range . 
B{B X,7, E m range) = — — — . (3) 

' IX JNg^X eKN460 

Figure |1] shows the efficiency-corrected spectrum using the KN460 model. We have taken a 
conservative approach of assigning to each bin a model-dependence uncertainty corresponding to 
the maximum efficiency deviation, among KN455 to KN480 models and the 12 BBU models, from 
the KN460 value. Note that these models span a range of mean energy up to 2a away from our 
mesured value. 

To make the results directly comparable to theoretical predictions, we correct for the fact that 
the cut in photon energy is made on reconstructed E* in the T(4S) (cm.) frame rather than on 
true E^ in the B rest frame. The latter energy has been smeared by a Doppler shift and by the 
(asymmetric) calorimeter resolution. (Doppler smearing arises from the Lorentz transformation 
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Figure 4: Efficiency-corrected photon energy spectrum for the extracted signal, shown only for the 
originally-blinded range of reconstructed energy (note the range 2.7-2.9 GeV is not used to measure 
the branching fractions or moments). The small error bar is statistical only. The larger error bar 
also includes BB and other systematic uncertainties and a model-dependence uncertainty, all in 
quadrature. There are significant correlations among non-statistical uncertainties for different bins. 



from the B rest frame to the T(4S') frame, given our ignorance of the direction of the B.) We 
compute the correction factor Ucut for the KN and BBU models and find that it has minimal model- 
dependence. Table|l]shows PBFs with corrections applied, along with the statistical, systematic and 
model-dependent errors. For the corrected PBFs the latter includes two correlated contributions: 
the model-dependent efficiency uncertainty noted above, already applied to the measured PBFs, 
and the uncertainty on acut- 

We studied many sources of systematic uncertainty, and here note the more significant. The 
uncertainty on the BB background subtraction is shown in Table|31 and amounts to 5.5% for 2.0 to 
2.7 GeV. It comes mostly from the statistical uncertainties on the correction factors derived from 
the vr'^(r/) control sample. Other systematic effects total 6.4% in quadrature. Of this, 3.3% is the 
uncertainty on photon selection, dominated by a 2.5% uncertainty on photon efficiency (determined 
from vr'^s in r decays) and 2% for the photon isolation cut. It also includes allowance for uncertainties 
in photon energy scale and resolution, and in the photon lateral shape cut efficiency, derived mainly 
from data from the BABAR B — >■ K*^ analysis and photons from virtual Compton scattering. The 
efficiency of the event shape cuts was studied using a vr*^ control sample to compare distributions of 
the Fisher discriminant between data and simulation, resulting in an uncertainty of 3.0%. A small 
sensitivity to details of Xg fragmentation implies, for the adjustments determined in the BABAR 
semi-inclusive analysis, an additional uncertainty of only 1.4%. A 2.2% uncertainty is assigned for 
lepton identification, and 3.0% for the uncertainties on the semileptonic corrections. 
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E range 
GeV 


Measured PBF (xlO~^) 
for reco E* range, cm. frame 


"cut 




Corrected PBF (xlO"^) 
for true range, B frame 


1.9 to 2.7 

2.0 to 2.7 

2.1 to 2.7 

2.2 to 2.7 


3.64 ±0.29 ±0.33 ±0.22 
3.32 ±0.26 ±0.28 ±0.16 
2.84 ±0.23 ±0.23 ±0.11 
2.23 ±0.19 ±0.18 ±0.06 


1.05 ±0.02 
1.07 ±0.02 
1.09 ±0.02 
1.13 ±0.03 


0.96 ±0.016 
0.96 ±0.016 
0.96 ±0.016 
0.96 ±0.016 


3.67 ±0.29 ±0.34 ±0.29 

3.41 ±0.27 ±0.29 ±0.23 
2.97 ±0.24 ±0.25 ±0.17 

2.42 ±0.21 ±0.20 ±0.13 



Table 4: Partial branching fractions (PBFs, preliminary) corrected to the B rest frame for four 
different reconstructed- £"* ranges. The uncertainties shown for the PBFs are statistical, systematic 
and model-dependent. The "Measured" values are in terms of E* in the cm. frame. The correction 
factor acut converts the measured branching fraction for an E* cut to a result in the B rest frame 
with an equal E^ cut, while the factor a^/^ removes the (4.0 ± 1.6)% B — > X^^'j component. For 
the last column, the error on acut is added linearly to the fractional model-dependence error, while 
the error on a^/s is added in quadrature to the systematic error. 

5.2 Truncated Moments 

As noted above, the truncated moments of the photon energy spectrum are of interest because 
they can be related to parameters of heavy quark theory. The measured first moment (E*) and 
second moment {E*"^) — {E*)"^ are computed by averaging over the efficiency-corrected photon 
energy spectrum in 100 MeV bins. We do this for four different energy ranges, extending from 1.9, 
2.0, 2.1 or 2.2 GeV to 2.7 GeV. The bin-by-bin efficiencies derived from the KN460 model are used 
to correct the yields. 

The observed moments are derived from the measured E* spectrum but the moments computed 
in theoretical calculations involve the true E^. Thus we need to correct the former for Doppler 
and calorimeter-resolution smearing. For a spectrum without energy cuts, the total shift of the 
first (second) moment is the sum of the first (second) moments of the smearing functions, about 
±0.005 GeV (±0.007 GeV^) for Doppler smearing plus -0.040 GeV (±0.014 GeV^) for the calorime- 
ter. The latter corrections are reduced when the spectrum is truncated by a minimum energy cut, 
e.g., to -0.027 GeV(-0.007GeV2) for a typical BBU model (m;, = 4.6 GeV, ^ o.45GeV2) with 
a cut of 2.0 GeV. An additional small shift accounts for cutting on the reconstructed E* in the 
T[AS) frame instead of on the true E^ in the B rest frame. For a given signal model, Monte Carlo 
simulation provides the total shift Atotai needed to obtain the first moment (-E^) from the measured 
first moment {E*) for the same numerical cut value; likewise for the second moment. 

TablesEland Elshow the observed and corrected first and second moment values for four different 
minimum photon energy cuts. Using 50 MeV bins instead of 100 MeV bins changes the first moment 
by less than 0.003 GeV (0.001 GeV) and the second moment by less than 0.001 GeV^ (0.0004 GeV^) 
for the 1.9 GeV cut (other cuts). The dominant systematic uncertainty arises from the BB back- 
ground subtraction. In addition the uncertainty in the photon energy scale and the modeling of 
photon energy resolution cause smaller systematic errors in the first and second moments, respec- 
tively. The model-dependent uncertainty on the first moment is based on redoing the computation 
using the efficiency for each KN or BBU model that matches the nominal first moment within 
its uncertainties. The second moment is sensitive to the modeling of the low-energy tail of the 
photon energy resolution but not to the efficiency model or to the uncertainty in the photon en- 
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E range 
GreV 


(E*) ± (stat) ± (sys) ± {model) 
GeV (with cut on E^) 


GeV 


(E^) ± (stat) ± (sys) ± (model) 
GeV (with cut on E^j 


1.9 to 2.7 

2.0 to 2.7 

2.1 to 2.7 

2.2 to 2.7 


2.270 ± 0.025 ± 0.017 ± 0.005 
2.304 ± 0.016 ± 0.010 ± 0.005 
2.350 ± 0.014 ± 0.007 ± 0.003 
2.412 ± 0.012 ± 0.005 ± 0.002 


0.018 ±0.007 
0.012 ±0.007 
0.005 ± 0.007 
-0.005 ±0.006 


2.288 ± 0.025 ± 0.017 ± 0.012 
2.316 ± 0.016 ± 0.010 ± 0.012 
2.355 ± 0.014 ± 0.007 ± 0.010 
2.407 ± 0.012 ± 0.005 ± 0.008 



Table 5: Measured first moments (prehminary) for various energy ranges, then corrected for resolu- 
tion smearing and for cutting on the B rest frame E^. The systematic errors have two sources: BB 
background uncertainties, with their full covariance matrix taken into account, and the uncertainty 
in the photon energy scale. The model-dependence error on {E*) is correlated with that on Atotah 
so these are combined linearly for the last column. The {E^) values in that column can be compared 
to theoretical predictions. Note that ~ 4% of i? — > X^j is included, but its first moments are not 
expected to differ significantly from those for B Xsj. 



E range 
GeV 


(E*^) -{E*f± Stat ± sys 
GeV^ (with cut on E*) 


^total 

GeV^ 


(E'^) — {Ej)'"^ ± stat ± sys ± model 
GeV^ (with cut on E^) 


1.9 to 2.7 

2.0 to 2.7 

2.1 to 2.7 

2.2 to 2.7 


0.0403 ± 0.0049 ± 0.0022 
0.0326 ± 0.0026 ± 0.0009 
0.0246 ± 0.0019 ± 0.0005 
0.0161 ± 0.0014 ± 0.0003 


-0.0075 ± 0.0025 
-0.0060 ± 0.0020 
-0.0055 ±0.0015 
-0.0045 ± 0.0005 


0.0328 ± 0.0049 ± 0.0023 ± 0.0025 
0.0266 ± 0.0026 ± 0.0010 ± 0.0020 
0.0191 ± 0.0019 ± 0.0006 ± 0.0015 
0.0116 ± 0.0014 ± 0.0004 ± 0.0005 



Table 6: Measured second moments (preliminary) for various energy ranges, then corrected for 
resolution smearing and for cutting on the B rest frame E^. The errors listed with Atotai represent 
model-dependence. The systematic errors have two sources: BB background uncertainties, with 
their full covariance matrix taken into account, and (after the correction) the uncertainty in the 
modeling of the photon energy resolution. The final corrected value for an E^/ cut can be compared 
to theoretical predictions. 

ergy scale. A systematic uncertainty for the photon energy resolution is determined by comparing 
the asymmetric photon resolution function measured using a sample of events from virtual Comp- 
ton scattering (in which the true photon energy can be determined by a kinematic constraint) 
with the corresponding function for simulated events. The uncertainty varies from 0.0005 GeV^ 
for E* > 1.9 GeV to 0.0002 GeV^ for E* > 2.2 GeV and is added in quadrature to the systematic 
uncertainty arising from the BB background systematic. We compute the corrections applied to 
the moments for the full range of KN and BBU models considered, in order to better represent 
possible model-dependence. The tabulated Atotai is the average value of the maximum and mini- 
mum corrections for these models, and its uncertainty is taken to be half the value of the difference 
between the maximum and minimum. The model-dependent uncertainties on {E*) and Atotai are 
correlated, and hence added linearly. 

Finally Tables [7| through |^ show the full correlation matrices between the eight measured 
moments for statistical, systematic, and model dependence errors, respectively. Statistical and 
systematic covariance matrices are computed using standard error propagation, starting from the 
measured spectrum in energy bins. For systematics, this incorporates the bin-to-bin covariance 
matrix that arises from the BB subtraction, with the fully-correlated covariance matrices arising 
from the energy-scale uncertainty (first moments only) and calorimeter resolution (second moments 
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only) added at the end. 

The model correlations take into account the variation of both the efficiency and the I^totai 
corrections among different models. For each of the models (indexed h) and each of the eight 
measured first and second moments Mj=i^8 we define A^g^y^^j(Mj) as the actual correction for that 
model. We also define Agjj(Mj) as the computed value of the first or second moment using a 
particular model minus its value using efficiencies from the nominal KN460 model. Letting 

we compute the covariance of A^^^^i{Mj) between each pair of moments by averaging over the mod- 
els, and thence derive the correlation matrix. Note that this approach, which treats the predictions 
of the different models as a statistical distribution, would yield smaller model-dependence errors 
than the ones we quote in Tables [5] and IHl but we use this result solely for the correlation matrix, 
not for the errors themselves. 







M M 


M 


M 


V 


V 


V 


V 


Quantity 


Min. 


1.9 2.0 


2.1 


2.2 


1.9 


2.0 


2.1 


2.2 


Mean 


1.9 GeV 


1.0000 0.5172 


0.3548 


0.2265 


-0.6110 


0.0077 


0.0821 


0.1008 


Mean 


2.0 GeV 


1.0000 


0.6838 


0.4326 


0.2285 


-0.0008 


0.1375 


0.1645 


Mean 


2.1 GeV 




1.0000 


0.6260 


0.4220 


0.5660 


0.1650 


0.1884 


Mean 


2.2 GeV 






1.0000 


0.4528 


0.7568 


0.7383 


0.2113 


Variance 


1.9 GeV 








1.0000 


0.4626 


0.3486 


0.2520 


Variance 


2.0 GeV 










1.0000 


0.6966 


0.4887 


Variance 


2.1 GeV 












1.0000 


0.6709 


Variance 


2.2 GeV 














1.0000 



Table 7: Correlation matrix (preliminary) for statistical errors on the truncated first and second 
moments (four each, for different minimum cuts). This matrix is symmetric, so redundant 
entries are not shown. Column labels are shorthand for the identical row labels. 







M M 


M 


M 


V 


V 


V 


V 


Quantity 


Min. E-y 


1.9 2.0 


2.1 


2.2 


1.9 


2.0 


2.1 


2.2 


Mean 


1.9 GeV 


1.0000 0.7841 


0.6076 


0.4622 


-0.7552 


-0.3862 


-0.1723 


-0.0942 


Mean 


2.0 GeV 


1.0000 


0.8007 


0.6311 


-0.2515 


-0.4753 


-0.2298 


-0.1451 


Mean 


2.1 GeV 




1.0000 


0.7886 


-0.0320 


0.0202 


-0.2903 


-0.2115 


Mean 


2.2 GeV 






1.0000 


0.0489 


0.1850 


0.2061 


-0.1001 


Variance 


1.9 GeV 








1.0000 


0.4473 


0.2219 


0.1436 


Variance 


2.0 GeV 










1.0000 


0.5003 


0.3342 


Variance 


2.1 GeV 












1.0000 


0.6592 


Variance 


2.2 GeV 














1.0000 



Table 8: Correlation matrix (preliminary) for systematic errors on the truncated first and second 
moments (four each, for different minimum E^ cuts). This matrix is symmetric, so redundant entries 
are not shown. Column labels are shorthand for the identical row labels. Systematic uncertainties 
included are those on the BB subtraction, the energy scale (first moments only) and calorimeter 
resolution (second moments only). 
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M M 


M 


M 


V 


V 


V 


V 


Quantity 


Min. 


1.9 2.0 


2.1 


2.2 


1.9 


2.0 


2.1 


2.2 


Mean 


1.9 GeV 


1.0000 0.9267 


0.9486 


0.8252 


-0.9057 


-0.9223 


-0.9234 


-0.7983 


Mean 


2.0 GeV 


1.0000 


0.9576 


0.8669 


-0.7452 


-0.8183 


-0.8087 


-0.8960 


Mean 


2.1 GeV 




1.0000 


0.9540 


-0.7390 


-0.7889 


-0.7843 


-0.7752 


Mean 


2.2 GeV 






1.0000 


-0.5378 


-0.5857 


-0.5788 


-0.6023 


Variance 


1.9 GeV 








1.0000 


0.9650 


0.9824 


0.6983 


Variance 


2.0 GeV 










1.0000 


0.9810 


0.8035 


Variance 


2.1 GeV 












1.0000 


0.8057 


Variance 


2.2 GeV 














1.0000 



Table 9: Correlation matrix (preliminary) for model-dependent errors on the truncated first and 
second moments (four each, for different minimum cuts). This matrix is symmetric, so redundant 
entries are not shown. Column labels are shorthand for the identical row labels. Uncertainties are 
included for the model-dependence of the signal efficiency and of the correction for resolution, 
Lorentz boost and cutoff. Sixteen models are averaged, as described in the text. 



The correlation matrices allow a fit of the moments to any parameterized theoretical calculation. 
However we note that of the eight moments only five are independent. One way to see this is by 
considering that there are two moments for the highest energy cut (E* > 2.2 GeV) and all the other 
moments can be derived by adding the three lower energy bins. Thus in fitting one must choose 
no more than five of the moments and check that the associated covariance matrices are positive 
definite. 

5.3 CP Asymmetry 

Within the Standard Model the combined direct CP asymmetry for 6 — )• 57 and b ^ dj decays, 

_ r(6 ^ S7 + 6 ^ dj) - r(b ^s^ + b^dj) . . 

~ r(6 ^ S7 + 6 ^ dj) + T(b^s-f + b^ dj) ' ^ ' 

is zero to a very good approximation New physics models with additional flavour violation 
can significantly enhance Acp to a few percent Furthermore, measurements of the combined 
Acp in 6 ^ 57 and 6 — > ^7 processes complement those of Acp in only b ^ s'j [201 El to 
constrain new physics models. 

The use of lepton tagging on the non-signal B for continuum suppression allows a straightfor- 
ward determination of Acp- The lepton's positive (negative) charge tags the signal side to contain 
a b (b) quark. Therefore 

- 1 
= iV+ + Ar-i_2.. ' 

where N^C) are the positive (negative) tagged signal yields and 1/(1 — 2uj) is the dilution factor 
due to the mistag rate, uj. 

The systematic uncertainty on the BB background subtraction only contributes a multiplicative 
uncertainty to the measured Acp- Therefore, the statistical uncertainty dominates and we find the 
optimal E* range for the Acp measurement to be 2.2 < E* < 2.7 GeV. The tagged signal yields are 

= 349 lb 48 and N~ = 409 it 45 where the uncertainties are statistical. We use control samples 
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of e^e^ — > X7 and B XTr^,rj to estimate the bias on Acp from any charge asymmetry in the 
detector or in the BB background. The bias is —0.005 it 0.013, which is apphed to the measured 
asymmetry. The uncertainty on the measurement of the bias is taken as an additive systematic 
uncertainty on the measured asymmetry. The measured asymmetry, uncorrected for dilution, is 
—0.084 zb 0.088 ± 0.013, where the first uncertainty is statistical and the second systematic. 

The mistag rate is (11.9±0.4)%. It is dominated by mixing in B^B^ events, which when averaged 
over B^B^ and B^B^ events is equal to half the time-averaged mixing rate = (9.3±0.2)% |22| . 
The remaining (2.6 it 0.3)% is estimated from simulation; these mistags are due to leptons from D 
decay, vr^ faking 7 conversions, vr'^ Dalitz decay, and charmonium decay. The error is dominated 
by uncertainties in the modeling of D decay and fakes. The dilution factor of 1/(1— 2ca;) = 1.31ib0.01 
scales the measured Aqp and its additive uncertainties to give 

Acp = -0.110 ± 0.115(stat.) ± 0.017(syst.) . (6) 

The systematic uncertainty includes a 5.4% multiplicative uncertainty from the BB background 
subtraction and 1.0% multiplicative uncertainty from the dilution factor. We have estimated that 
any model dependent uncertainty, due to differences in the B — >■ Xg'y and B X^j spectra, is 
much smaller than the systematic uncertainty. 

6 Conclusions 

We have performed a lepton-tagged fully-inclusive measurement oi B ^ Xg'j decays. We present a 
reconstructed photon energy spectrum in the T{4S) frame, and partial branching fractions above 
minimum reconstructed energies of 1.9, 2.0, 2.1 and 2.2 GeV (Table|lJ). We also present the partial 
branching fractions, the first and second moments of the true photon energy distribution in the B 
rest frame and the correlations between the first and second moments, above the same minimum 
photon energies (tables El El El 13 El and 121). Finally, we have used the charge of the tagging leptons 
to obtain the CP asymmetry of the sum of i? — > Xsj plus B X^^l with reconstructed photon 
energy above 2.2 GeV. 
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We present preliminary results from a lepton-tagged fully-inclusive mea- 
surement of B ^ Xg^ decays, where Xs is any strange hadronic state. 
Results are based on a BaBar data set of 88.5 million BB pairs at the 
T(45') resonance. We present a reconstructed photon energy spectrum in 
the T(45') frame, and partial branching fractions above minimum recon- 
structed photon energies of 1.9, 2.0, 2.1 and 2.2 GeV. We then convert these 
to measurements of partial branching fractions and truncated first and sec- 
ond moments of the true photon energy distribution in the B rest frame, 
above the same minimum photon energy values. The full correlation matri- 
ces between the first and second moments are included to allow fitting to 
any parameterized theoretical calculation. We also measure the direct CP 
asymmetry Acp{B — > X^j^dl) (based on the charge of the tagging lepton) 
above a reconstructed photon energy of 2.2 GeV. 



